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МАШИНОСТРОЕНИЕ 


УДК 681.51 
В.Л.ЗАКОВОРОТНЫЙ *, НГУЕН ДОНГ АНЬ””, ФАМ ДИНЬ ТУНГ* 


УСТОЙЧИВОСТЬ ЭВОЛЮЦИОННОЙ ТРАЕКТОРИИ 
МЕХАНИЧЕСКОЙ СИСТЕМЫ, . 
ВЗАИМОДЕЙСТВУЮЩЕЙ С ТРИБОСРЕДОЙ 


Рассматриваются условия потери устойчивости движения механической системы, 
взаимодействующей со средой, формируемой в трибосопряжении. Показаны усло- 
вия формирования и роль диссипативных, гироскопических, потенциальных и цир- 
куляционных сил, формируемых в узле трения естественным образом и изменяю- 
щихся в ходе эволюционных преобразований. 

Ключевые слова: эволюция, трибосреда, устойчивость, стационарная траекто- 
рия. 


Введение. Для обеспечения движения элементов машин и механизмов во 
многих случаях необходимо решать вопросы устойчивости движения. При 
этом, как правило, некоторый элемент машины взаимодействует с узлом 
трения. В свою очередь, в зоне, находящейся между контактируемыми по- 
верхностями, формируется диссипативная среда, образующая динамиче- 
скую связь [1-3]. Основные динамические свойства такой системы можно 
раскрыть на основе использования «базовой» динамической модели трибо- 
системы [4], рассматривающей смещение индентора, подвешенного с по- 
мощью упругодиссипативных подвесок к абсолютно жёсткому основанию. 
Индентор в базовой модели представляется в виде сосредоточенной массы. 
Поэтому изгибные колебания его поверхности не учитываются. В свою оче- 
редь, образец также обладает большой массой и считается недеформируе- 
мым (рис.1). Он движется со скоростью Г = соп5Е , задаваемой от внешне- 
го источника. Тогда уравнение движения такой системы можно предста- 
вить в виде 





„4х, „ах ах 
+ ры СХ = ВХ, ,Р)+0, (1) 
Ей а 
ге Х={Х,,Х,} - вектор состояния системы, координаты которого 


показывают текущие значения положения поверхности индентора в 
трибосреде (формируемой в процессе относительного скольжения 
переходной области между контактируемыми поверхностями); 


[Хх Хх [Хх 
в(х,й ‚Г,р)= (вой: БЕ,“ ‚Г.Р - вектор-фун- 


кция, раскрывающая зависимость сил контактного взаимодействия от 
координат состояния, её свойства определяются формируемой в 
процессе трения переходной областью между индентором и 
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образцом, названной трибосредой; р= {р,р., Р,} - параметры, 
конкретизирующие нелинейные —вектор-функции; ее. - 
положительно определённые симметричные матрицы размерностью 
282 соответственно обобщенных масс, — коэффициентов 
демпфирования и жесткости. В рассматриваемой модели 
дополнительно 7 является диагональной; И- вектор внешних 
силовых воздействий. В статье рассматривается случай, когда 
И = соп5Ё и Г-= соп5Ё. Здесь и далее символ +$--}” означает 
операцию транспонирования. 














Основание 


крепления 
индентора 


Трибосреда ы 






__ Е= {Е (ХЕХ) 


Рис.1. Схема «базовой» динамической модели узла трения 


Трибосреда и, следовательно, параметры Р формируются в 
результате работы сил контактного взаимодействия, и свойства 
трибосреды определяются текущим значением мощности необратимых 
преобразований в контактной области с учётом её предыстории изменения. 
Следовательно, свойства трибосреды зависят от траектории мощности 
необратимых преобразований в контактной области по совершаемой 
работе. Поэтому параметры Р можно представить в виде интегрального 


оператора 
А 


А 
В (0 (2) а) 9 <) 
Р(А)= ро+в; Ача; № <-АМс)ча; м (4-9)МС 4 , (2) 
0 0 
1 2 
где и? (с - А), и” ($ - 4) - ядра интегральных операторов; хо - 
2 Е) 
исходное значение параметра; © ‘”, а, @аС’- коэффициенты; 
А, М№ - работа и мощность необратимых преобразований в зоне 
контакта, то есть №(.4) - фазовая траектория в плоскости «работа — 
мощность». 

Ядра интегральных операторов моделируют две противоположные 
тенденции в изменениях параметров системы. Первый интегральный опе- 
ратор характеризует процесс стабилизации свойств контакта и выхода ди- 
намической системы на некоторый стационарный режим. Это установление 
равновесной шероховатости, стационарного трибохимического обмена 


между поверхностями контактируемых тел и трибосредой, установление 
стационарного гидродинамического процесса в контактной области в слу- 
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чае граничного трения и пр. Этот процесс характеризуется, как правило, 
уменьшением переходной области между контактируемыми поверхностями. 
Второй интегральный оператор характеризует деградацию контактной об- 
ласти и в итоге — установление режима катастрофического изнашивания 
поверхностей, влияющего на параметры Р. Он характеризуется непропор- 
ционально быстрым увеличением переходной области -— трибосреды. 

В статье рассматривается случай, для которого справедливо 


Г АХ, / 4, то есть изменение знака текущей скорости относительного 
скольжения отсутствует. Наиболее важные свойства вектор функций 
Е(Х, АХ /4ЕТ,р) для этого случая заключаются в следующем [1-3]: 

- ВО, аХ/аЬУТ,р)= Е(Х,Г,р), то есть нормальная состав- 


ляющая силы контактного взаимодействия зависит от сближения поверхно- 
стей, и она не зависит от стационарного смещения индентора относитель- 
но образца в тангенциальном направлении; 

- ЕСсаХ/а,Т,р)= Е[У (АХ, Ш+И),Р], У= Х(-Г), о 
есть изменение тангенциальной силы запаздывает по отношению к измене- 
нию её нормальной составляющей. Кроме этого тангенциальная сила зави- 
сит от текущего значения скорости относительного скольжения и возмож- 
ны её значения, при которых увеличению скорости соответствует уменьше- 
ние силы №. 


Параметры Р фактически есть коэффициенты аппроксимации 
рассматриваемых нелинейных функций. Они обладают свойством эволюци- 
онной изменчивости. Однако в пределах импульсных реакций системы их 
можно считать постоянными, а эволюционную систему при заданной рабо- 
те — замороженной. В связи с этим имеется возможность ввести в рассмот- 
рение понятие стационарной эволюционной траектории Х” (4), которая 
удовлетворяет следующим требованиям: 


сХ' = Е(Х',0,Г,р)+И; (3) 


А А 
РСА) = риа АЗЯТУ с - ДЕС) ча м” (4-6)Е.6 46 ‚ (4) 
0 0 
где р,(А) - эволюционная траектория Г - го параметра. Следовательно, в 
параметрическом пространстве р“ мы имеем медленно изменяющуюся 
траекторию р(А)Р *, которой соответствует медленно меняющаяся эво- 
люционная траектория Х” (А). В ходе эволюции система (1) может поте- 


рять устойчивость. Тогда в окрестности эволюционной траектории Х” (А) 


формируются многообразия (орбитально асимптотически устойчивые пре- 
дельные циклы, инвариантные торы, хаотические аттракторы). Однако для 
их анализа необходимо раскрывать нелинейные вектор-функции 


ах 
ЕСХ и Р), что будет сделано в следующих наших публикациях. 
Здесь важно подчеркнуть, что с помощью уравнений (3) и (4) можно вы- 


числить траектории Х“(А) и р(А)Р *, а затем уже анализировать их 
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устойчивость. Причём, траектории Х“’(А) и р(А)Р * всегда существу- 
ют, так как уравнение (3) характеризует квазистатику. 

Для изучения же условий потери устойчивости Х” (А) необходимо 
рассмотреть уравнение в вариациях относительно этой траектории. Тогда 
































для вариаций | ^® | то есть | х<«р—х`-хе> | имеем систему в 
вариациях относительно |х |: 
ах. м ах 
т + сх= (Хх, —), (5) 
Аг @ [ИА 
где 
4х 


о: ЕСА + х)- АХОЫЕКХ, + у), “- ВХ) = Ф.О, 


у= х (1-1) - нелинейные функции в вариациях относительно точки 


равновесия, причем © (Х’,0,0) 0. Устойчивость траектории Х’(А) 
определяется свойствами автономной системы 

















Я?х Ах 
в +Й; ре , (6) 
т 
аа (си+ с) Сл С, й й 1 й, , а 90 а , 
= т) - (Т) (г) м то 1 1 
(ета) 652 ' й,.-й. ) (й,>- р, | } 


с = -90,/9м; Йо =190,/9х; #0 =90,/9х,. В (6) принято во 
внимание, что 90,/9х 0, 960,/9х 0 и рассматривается скоростной 


диапазон, в котором 90 ,/9х, О. Все указанные параметры зависят от 


внешних сил С/ ‚ скорости относительного скольжения Г’ и внутренних 
факторов, определяемых работой А. Таким образом, они играют роль 
управляющих параметров, перестраивающих динамические свойства 
системы. В рассматриваемой системе роль управления играет работа сил 
контактного взаимодействия, изменяющая следующие параметры: 


Т = 
-В = Е, - динамическую контактную жесткость трибосопряжения; 
_ (Г) о 
-Р› = С,› - параметр, зависящии от динамического коэффициента 


трения А=90,/9д0,, то есть ты с’. Заметим, что динамический 

коэффициент трения может существенно отличаться от традиционно 
Т 

рассматриваемого коэффициента трения К“? = Е. (Х\)/Е(Х,), так как 


он связывает нормальные и тангенциальные составляющие сил в 
вариациях относительно их значения в точке равновесия, 


-Рр.=Т - величину запаздывающего аргумента в преобразовании 
вариаций нормальных составляющих сил в тангенциальные; 
Т 
=Ра= о > коэффициент вязкого сопротивления в тангенциальном 


направлении. Он может иметь и отрицательное значение в скоростном 
диапазоне, где увеличению скорости относительного скольжения 


соответствует уменьшение силы трения Я . 
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Прежде всего остановимся на структуре динамических сил, 
влияющих на устойчивость траектории Х” (А). В свою очередь, структура 
динамических сил зависит от эволюционной траектории параметров 
Р(А)Р ‘. 

Структура динамических сил. Её влияние на устойчивость 
эволюционной траектории. Исходные динамические подсистемы без 
трибосреды, задаваемые матрицами р Е (левые части в (5)), 
имеют постоянные параметры. Эти матрицы, как уже отмечено, являются 
положительно определенными и симметричными, так как динамическая 
структура машины, как правило, является неизменной, а силовые функции 
обладают потенциальными свойствами. Кроме этого исходная система без 
трения имеет асимптотически устойчивую точку равновесия. Напомним, 
что диссипативная функция при моделировании динамики машин вводится 
в форме Релея. Однако при переходе к (6) суммарные матрицы диссипации 

































































‚ | И жёсткости | ‹ | уже не обладают симметричными свойствами. 
Поэтому с учётом (6) их можно представить в виде суммы [5-8] 
еее к р И; = И + и р (7) 
(т) (т) т 
где с® = . (са + с ) (с, 2+ 0, 5 >) (с) _ т (> Е Оо 
(т) ит т (Т 
(5+ 0,565”) в (и ›- 0,51 (2- Но 
- симметричные части матриц жёсткости и диссипации, отвечающие 
: ® 0 - 0,5612 
за потенциальные свойства системы; с;’= (г) ; 
0,5, 0 
2 0 0,58 (к) 
"= - кососимметричные матрицы. Матрицы с,>, 


(т) 
>05 0 
как известно, формируют непотенциальные (циркуляционные) силы. 
Они могут привести к потере устойчивости точки равновесия. 


Матрицы И) формируют гироскопические силы, которые могут 


стабилизировать точку равновесия. 

Мы приходим к важному заключению: 

1. При взаимодействии механической системы с трибосредой 
естественным образом формируются гироскопические и циркуляционные 
силы, влияющие на условия асимптотической устойчивости системы. 
Подчеркнём, что все имеющиеся экспериментальные данные по изучению 
траекторий движения контактируемых поверхностей в узлах трения 
показывают, что при потере устойчивости траектории движения индентора 
всегда остаются круговыми независимо от параметров подвески индентора. 
В этом подтверждение роли циркуляционных сил, то есть сил, 
ортогональных направлению упругих деформаций. 











Ра Т 
2. Симметричная часть матрицы | * | за счёт коэффициента м, 


определяющего отрицательное приращение тангенциальных сил трения 
при увеличении скорости относительного скольжения, может 
преобразоваться из положительно определённой в отрицательно 
определённую. Таким образом, в системе дополнительно могут 
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формироваться ускоряющие силы, приводящие к потере устойчивости 
точки равновесия. 

Для определения влияния структуры сил на устойчивость 
представим (6) в виде 


ах Хх Е к) ах : 
> ое И + сх=0. (8) 
Е @ Гоа 
Рассмотрим несколько случаев. 


Первый случай. Рассмотрим влияние асимметрии позиционных сил 
Т 
в предположении, что т > би #5. - 0. Тогда #® - 0, 1? = и =йи 
вместо (8) имеем 





+7 


ах Ах р р 
НИ с х+ сх= 0. (9) 
Ф И 
Мы видим, что позиционная Связь, формируемая процессом трения, 


вызывает изменение потенциальных свойств динамической системы и 
приводит к образованию циркуляционных сил (характеризуются матрицей 


2” ), которые могут влиять на устойчивость точки равновесия. Из (9) 


получаем характеристический полином системы 

А; (р)=А(р)+ с [тр® + (#,)р+ (се, Л С ( 6: + 21) ‚ (10) 
где А(р) - характеристический полином исходной системы без трения. 

Полином А (р) соответствует асимптотически устойчивой системе. 

Анализ (8) позволяет сформулировать первое свойство влияния 
асимметрии позиционных сил, формируемых в трибосопряжении, на 
устойчивость точки равновесия системы. Если исходная динамическая 
структура подвески индентора является ортогональной, то 
циркуляционные силы не могут привести к потере устойчивости 
равновесия. Здесь и Далее под ортогональной будем понимать 
динамическую структуру, для которой все недиагональные элемента 
матриц Сий равны нулю, то есть в исходном СОСТОЯНИИ 
координаты Хх ЯВЛЯЮТСЯ нормальными. Действительно, при 


условии с.=й.=0О полином (10) преобразуется в полином 
= 2 (Г) 2 

А;(р)= (тр + й.р+ Сл + С )(тр` + й.р+ С›2), корни которого 

имеют отрицательные вещественные части. Очевидно, что это 

утверждение можно распространить и на случай, когда механическая часть 

системы без трения представлена в виде № -мерной динамической 


структуры, обладающей ортогональными динамическими свойствами. 


Т о 
Обозначим 0:55) = И и рассмотрим характеристический полином 





т 


системы, приведённой к виду (9), при й“® = 0, Й= А: 
А, (р)=А(ру+и’, (11) 


-. 2 (с) 2 (с) (<) \2 
где А (р)= (тр + йр+ С {тр` + й,р+ с›2)- (И ›р+ с) 
характеристический полином системы (11) без связей, формирующих 
циркуляционные силы, то есть обладающей потенциальными свойствами. 


Пусть матрица со является положительно определённой. Матрица Й 
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является положительно определённой. Поэтому полином А (р) имеет все 


корни с отрицательными вещественными частями, то есть соответствует 
асимптотически устойчивой системе. Для выяснения ВЛИЯНИЯ 


циркуляционных сил, учитываемых д ^, удобно воспользоваться частотным 


критерием устойчивости Михайлова (рис.2). Очевидно, что и” смещает 
годограф Михайлова вдоль вещественной оси во всём частотном диапазоне 
на величину и” и может преобразовать асимптотически устойчивую 


систему в неустойчивую. Кроме этого при увеличении Й значение 


я Т 
определителя, составленного из матрицы т по мере увеличения бо 


уменьшается. Принципиально при неограниченном увеличении К матрица 


в может стать отрицательно определённой. 





Рис.2. Схема преобразования годографа Михайлова 
за счёт циркуляционных сил, формируемых в трибосопряжении 


Мы видим, что неконсервативные (циркуляционные) силы, 
естественным образом формируемые в любой динамической системе 
трения, могут приводить к потере устойчивости точки равновесия системы. 
Причём потеря устойчивости зависит от параметров динамической 
подсистемы индентора, характеристик трибосопряжения, внешних условий 


и от траектории р(А)Р ^. Подчеркнём, что положение равновесия 


определяет точку, в окрестности которой определяется градиент функции 
сближения и динамический коэффициент трения. 

Второй случай. Рассмотрим влияние на устойчивость равновесия 
запаздывающего аргумента т, изменяющего симметричную 


(потенциальную) матрицу — диссипации Я и формирующего 
гироскопические силы за счёт кососимметричной части й;^®. Как и в 


первом случае, если матрица диссипации исходной системы является 
диагональной (подвеска индентора — обладает — ортогональными 
динамическими свойствами), то запаздывающий аргумент никак не влияет 
на устойчивость точки равновесия системы. В общем же случае структура 
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К ь 
матриц д и р ) показывает, что запаздывающии а мент Т 
х х 


усиливает положительную определённость матрицы р и дополнительно 


приводит к формированию гироскопических сил, также увеличивающих 
запас устойчивости системы. Для доказательства этого утверждения 
рассмотрим характеристический полином системы в предположении, что 


1’ = 0. Тогда имеем характеристический полином системы АД”. (р): 
А (РЕ А, (р)+ ТАС) р(й..Р + с,1), (12) 
где А; (р) - соответствует (10). 

Пусть полином А; (р) соответствует системе, устойчивой по 
Ляпунову. Тогда годограф ЛД; (74) проходит через начало координат 
(рис.З‚а). Каждая точка годографа Л’, (10) смещается влево на 
— Ас’ й..@” и вверх на величину 1Ас\’с,.@7, то есть будет 
смещаться во внешнюю сторону по отношению к годографу А; (р), 


соответствующему устойчивой по Ляпунову системе. Следовательно, 
годограф А’, (10) будет отвечать асимптотически устойчивой системе. 


Более того, если годограф Л’, (70) соответствует неустойчивой системе, 


то по мере увеличения Т система может приобрести асимптотическую 
устойчивость (рис.3,6). Подчеркнём, что механизм влияния Т на 
устойчивость равновесия двоякий. Во-первых, по мере увеличения 
запаздывающего аргумента усиливается положительная определённость 
потенциальной матрицы диссипации, во-вторых, при этом возрастает 
стабилизирующее влияние — гироскопических сил, — формируемых 


кососимметричными составляющими матрицы д . 


Сделаем вывод.  Гироскопические члены, — определяемые 
запаздыванием  тангенциальной составляющей силы — контактного 
взаимодействия по отношению к смещениям индентора относительно 
образца в нормальном направлении, стабилизируют точку равновесия. 





Рис.3. Влияние на годограф Михайлова запаздывающего аргумента 7 
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Третий случай. Рассмотрим влияние на устойчивость параметра 
Т К К у 
№7, который связан с кинетической характеристикой процесса трения 


[4]. Нетрудно показать, что к аналогичным эффектам приводит и запазды- 
вание нормальных сил по отношению к нормальным к поверхности смеще- 


ниям индентора по отношению к образцу. Анализ матриц Я; и с; пока- 


Г 
зывает, что параметр В влияет на симметричную составляющую п Е 


Т 
Более того, при увеличении йо ВОЗМОЖНО преобразование матрицы 


Ио в отрицательно определённую. В этом случае в системе трения могут 


формироваться ускоряющие силы, при всех условиях вызывающие потерю 
устойчивости точки равновесия системы. 

Для выяснения механизма потери устойчивости В этом случае 
вновь воспользуется изложенным выше приёмом. Для этого рассмотрим 


характеристический полином А! (р): 

А; (р)=А,(р)- у ртр” +7 Р+си+ И (13) 
и выясним, каким образом смещается годограф Михайлова за счёт члена, 
зависящего от во (рис.4), если полином Д'(р) соответствует системе, 
устойчивой по Ляпунову. В терминологии теории управления система, 
имеющая полином А’ (р), находится на границе устойчивости. 





Рис.4. Преобразование годографа Михайлова 
за счёт влияния коэффициента отрицательного динамического трения 
в тангенциальном направлении 


Подчеркнем, что формируемые в контакте гироскопические и 
циркуляционные силы в этом случае не могут стабилизировать равновесие 
системы, если оно потеряло устойчивость за счёт отрицательного 
динамического трения или за счёт запаздывания нормальных 
составляющих сил по отношению к их нормальным смещениям. 
Удивительным является то, что при увеличении контактной жёсткости 
системы в нормальном направлении, а также при возрастании жёсткости 
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подвески и вязкого сопротивления индентора в этом же направлении 
склонность К потере устойчивости равновесия системы возрастает. 
Влияние на устойчивость эволюционных параметров. Так как пара- 


метры трибосреды р,(/А),1= 1,2,3,4 обладают свойством эволюционной 


изменчивости, то при анализе устойчивости медленно изменяющейся эво- 
ЛЮЦИОННОЙ траектории в целом в пределах всего периода работы сил тре- 
ния естественно воспользоваться методом Ш -разбиения. Условие устойчи- 


вости в целом определяется требованием р(А) М (М - область 


устойчивости в пространстве варьируемых параметров). 

Проиллюстрируем анализ устойчивости эволюционных траекторий 
на примере, анализ которого имеет самостоятельное значение. Рассмотрим 
систему, параметры динамической характеристики процесса трения в кото- 
рой удовлетворяют следующим требованиям: 1 - 0, на индентор прило- 


жена внешняя нормальная сила Р, = соп5Ё . Кроме этого рассмотрим сум- 

марный коэффициент вязкого трения в тангенциальном направлении 
# Т о Т 

№о= Ия Жо ‚ который может в зависимости от й_) = р. принимать от- 


рицательное значение. Так как динамический коэффициент трения может 
существенно отличаться от традиционного коэффициента трения, введём в 


Т К о 
рассмотрение параметр ру = К?) , обладающий также свойством эволюци- 
онной изменчивости, то есть № = р.ЁН. В рассматриваемом примере 
Йсон51 . В связи с тем, что р. тоже эволюционный параметр, рассмот- 


рим совокупность областей устойчивости для различных значений р.. 


Экспериментальные данные, приведённые в работах [1-3], показы- 
вают, что традиционная система, функционирующая в режиме трения без 
смазки или в случае граничного трения, имеет две тенденции. 

Первая направлена на стабилизацию равновесия. При этом после 
приработки возрастает контактная жёсткость трибосопряжения в нормаль- 


ном направлении (параметр р, ); уменьшается динамический коэффициент 
трения, следовательно, параметр Р,›; уменьшается коэффициент трения 
(параметр Ру); практически остаётся неизменным величина запаздываю- 
щего аргумента (параметр рР.), так как он в основном зависит от средней 
скорости относительного скольжения (в нашем примере р) = 0); уменьша- 
ется коэффициент отрицательного вязкого трения в тангенциальном 


направлении (параметр р. ). 


Вторая тенденция направлена на деградацию контактной области. 
Она проявляется в противоположных тенденциях изменения указанных па- 
раметров. 


В связи с этим для определения эволюционных параметров р и 


Р›, изменения которых зависят от эволюции Ру, можно воспользоваться 
следующей системой: 
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0 1 


0 1 


} -А ИЕ 
р.(4д= рота А-а Ру окно ор. 
0 0 2 


2 





Первый интегральный оператор позволяет определить эволюцион- 
ную траекторию изменения тангенциальной составляющей силы трения 


Е, (А) ‚ остальные уравнения позволяют определить эволюционные траек- 


тории параметров р,(А), р.(А). Для определения области варьируемых 
параметров, при которых система асимптотически устойчива, естественно 
использовать метод ЛД -разбиений [9]. Полагая р. › = СОП5Т , из (6) получа- 
ем уравнения фигуративных линий: 

(- тб 2+ сз) р - а2рР> = - [4 0 “- 4 ›0 +0 


т (15) 
й, ›0 Р- й, 20 р. = -0 [@;0 2 Ы 


> 2. = ы _ : 
где Чи = С. 1С>›- (с,2) ; @, = Иа + а 2Й 5615 р 


— } 2 . ыы у . — 2 
@ = пси + с,,)+ Ай, - (И >); @,= т(йи+ №,); @.= т”. Система 
(15) позволяет построить фигуративные линии и выделить область устой- 
чивости системы. Задаваясь дискретными значениями Р, на отрезке 


(0. (2) м м и 
р. (Р.Р. ), получаем семейство областей устойчивости. Пример об- 


ласти устойчивости показан на рис.5 для системы, имеющей следующие 


1000кГ'/ мм; 1400кГ/ мм 


параметры подвески индентора. <= Е. 2250кГ/ мм 7 


0,1кГс/ мм; О13кГс/ мм 
0,13к/Гс/ мм; 0,2кГс/мм 


_ | 10? кГс? / мм; 0 
0; 10 кГс? / мм 


гральные операторы (14), приведены в таблице. 


| Параметры, входящие в инте- 


Параметры, входящие в интегральные операторы 


ее м 





[0 | 2.105 [кГм]* | 1,2-10°, сДкГм]р | 5-10, с/[кГм 
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20, м? 120, с/[кГм? 0,3, с/[кГм? 500 4000 


1 
400, скГм?] | 0,5, сДкГм?] 3200 











4 
Р.Г 
Рис. 5. Области О-разбиения, в которых приведена эволюционная траектория 
параметров р,(.А) и р.(А) при изменении А в пределах (0;1,8 10*кГм). Коор- 
динаты точек: А(р = 2 10°, р,=2 10); В(р, = 2,14 10°, р, = 3,06 10°); 
СС = 2,07 103, р» = 3,56 10%) 


Прежде всего охарактеризуем область устойчивости системы. Для 
р. = 0 (коэффициент вязкого трения со стороны трибосреды равен нулю) 
приведена двойная штриховка фигуративной линии и одинарные штрихов- 
ки особой линии при @ = 0. Таким образом, область устойчивости нахо- 


дится ниже приведённых линий для р. = 0. По мере увеличения р. об- 


ласть трансформируется, как показано на иллюстрации, и при р. 0,031 


вырождается полностью. В этом случае матрица диссипации с учётом влия- 
ния трибосреды становится отрицательно определённой. 
На область устойчивости варьируемых параметров наложена эво- 


люционная диаграмма р!(4) и р,(4А), которая при р. = 0 в точке «с» 
пересекает одинарную штриховку, а при р.= 0,01 - вначале двойную 
штриховку, а затем одинарную. Этим двум системам, отличающимся значе- 


нием р. = 0 и р. = 0,01, соответствует эволюционная диаграмма двух 


о 1 2. о 
корнеи ро и р . характеристического полинома замороженнои системы 
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(рис.6). Два других корня, имеющих значительные отрицательные веще- 
ственные части, на рис.6 не показаны. Таким образом, в ходе эволюции си- 


стемы при р. = 0 система в точке «с» приобретает свойство апериодиче- 
ской неустойчивости. При р. = 0,01 она в точке «Б» становится неустой- 
чивой, причём неустойчивость имеет колебательный характер. 





Рис.6. Эволюционные траектории корней характеристического полинома системы 


Для рассмотрения эволюции системы в области, в которой эволю- 
ционная траектория является неустойчивой, необходимо привлекать нели- 
нейные уравнения связи, формируемые трибосредой. В настоящей статье 
этот вопрос не рассматривается. Общие свойства рассматриваемой в каче- 
стве примера системы таковы: на первом этапе приработки системы трения 
запас устойчивости в системе возрастает. В пределах работы сил трения 


(0,15 10*;1,2 10*)кГм динамические свойства системы остаются пример- 


но неизменными. Затем система достаточно быстро теряет устойчивость 
эволюционной траектории, и в зависимости от параметров подвески инден- 
тора потеря устойчивости может быть апериодической или колебательной. 
Заключение. 1. При взаимодействии механической системы с трибо- 
средой получаемая динамическая система обладает принципиально иными 
динамическими свойствами по сравнению с исходной механической систе- 
мой без трения. При этом в контакте в зоне трения естественным образом 
формируются диссипативные, ускоряющие, гироскопические, потенциаль- 
ные и циркуляционные силы. Характерно, что диссипативные и потенци- 
альные силы зависят как от параметров связи, формируемой узлом трения, 
так и параметров подвески индентора. Что касается гироскопических 
и циркуляционных сил, то они не зависят от параметров подвески инденто- 
ра. 

Механизм изменения потенциальных и формирования циркуляцион- 
ных сил связан с тем, что позиционные смещения контактируемых поверх- 
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ностей в направлении скорости относительного скольжения не изменяют 
силовых реакций со стороны трибосопряжения. При определённом соотно- 
шении нормальных и тангенциальных составляющих сил контактного взаи- 
модействия, рассматриваемых в вариациях относительно стационарного 
состояния, динамическая система за счёт формирования циркуляционных 
сил может потерять устойчивость. Потеря устойчивости зависит также от 
параметров подвески индентора. В частности, если матрицы жёсткости и 
диссипации подвески индентора являются диагональными, то циркуляци- 
онные силы не влияют на устойчивость равновесия. В этом случае возмож- 
на потеря устойчивости тангенциальной формы колебаний за счёт отрица- 
тельного коэффициента трения в вариациях относительно стационарного 
состояния. Этот случай фактически вытекает из известных моделей Релея и 
Ван дер Поля, которые после некоторой модификации используются для 
объяснения формирования автоколебаний механических систем, взаимо- 
действующих с трением. 

Механизм изменения потенциальных диссипативных и формирова- 
ния гироскопических сил обусловлен влиянием запаздывания изменения 
тангенциальных составляющих сил по отношению к вариациям нормальных 
составляющих смещений индентора относительно образца. Это запаздыва- 
ние всегда направлено на стабилизацию точки равновесия системы. За 
счёт этого система, устойчивая по Ляпунову, всегда становится асимптоти- 
чески устойчивой. Более того, при увеличении запаздывания неустойчивую 
точку равновесия можно стабилизировать за счёт изменения потенциаль- 
ной матрицы диссипации и формирования гироскопических сил. Влияние 
этих сил на устойчивость равновесия зависит от параметров подвески ин- 
дентора. Причём, гироскопические силы не влияют на устойчивость равно- 
весия в том случае, если недиагональные элементы матриц жёсткости и 
диссипации подвески равны нулю. 

Механизм формирования ускоряющих сил обусловлен существова- 
нием кинетической характеристики процесса трения, то есть уменьшением 
тангенциальной силы трения по мере увеличения скорости относительного 
скольжения в определённом скоростном диапазоне. Ускоряющие силы мо- 
гут формироваться и за счёт запаздывания нормальных к контактируемой 
поверхности сил при вариациях смещения поверхности в этом же направ- 
лении. 

2. Параметры линеаризованной динамической связи, формируемой 
трибосопряжением, обладают свойством эволюционной изменчивости. 
Причём изменения параметров определяются траекторией мощности 
необратимых преобразований в узле трения по совершённой работе. 
Поэтому текущие значения параметров можно представить в виде инте- 
гральных операторов типа операторов Вольтера второго рода. Как прави- 
ло, изменения параметров имеют две тенденции. Первая направлена на 
стабилизацию медленно изменяющейся эволюционной траектории. В этом 
случае запас устойчивости в системе возрастает. Вторая — на деградацию 
свойств контакта. Эта тенденция вызывает потерю устойчивости. Подчерк- 
нем, что выполненный анализ относится к малым колебаниям относительно 
точки равновесия, поэтому при потере устойчивости в системе необходимо 
анализировать нелинейные связи, формируемые процессом в вариациях 
относительно точки равновесия, и образующиеся при этом многообразия в 
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пространстве состояния системы, которые являются естественными для 
рассматриваемой системы. 

3. Если задана априорная информация об изменениях внешних 
условий функционирования системы (функция изменения нормального дав- 
ления и скорости относительного скольжения), то предложенные матема- 
тические модели позволяют прогнозировать динамическую перестройку си- 
стемы. Кроме этого на основе наблюдения за траекториями корней харак- 
теристического полинома системы имеется возможность осуществлять ди- 
намический мониторинг состояния узла трения в процессе эксплуатации. 
Одним из направлений мониторинга является авторегрессионный спек- 
тральный анализ наблюдаемых вибрационных последовательностей. 

В статье излагаются результаты исследований, полученных при 
поддержке гранта 07-09-90000 РФФИ по Российско - Вьетнамскому научно- 
му сотрудничеству. 


Библиографический список 

1. Заковоротный В.Л. Нелинейная трибомеханика / В.Л. Заково- 
ротный. — Ростов н/Д: Издательский центр ДГТУ, 2000. - 293 с. 

2. Заковоротный В.Л. Динамика трибосистем. Самоорганизация, 
эволюция / В.Л. Заковоротный. — Ростов н/Д: Издательский 
центр ДГТУ, 2003. - 502 с. 

3. Заковоротный В.Л. Введение в динамику трибосистем / В.Л. За- 
коворотный. -— Ростов н/Д: ИнфоСервис, 2004. - 680 с. 

4. Крагельский И.В. Фрикционные автоколебания /И.В.Крагель- 
ский, Н.В.Гитис. — М.: Наука, 1987. 

5. Иеви С. иг ВазНКайоп уоп КгаКеп. Зсймем / САМепт, 
Н.Дедег. Ваигекипд, 84, № 48, 1966. 

6. Дебег Н. Мпеаг Еа$Нс За Ку / Н.Иедег. СИЧса! Апау$5 оЁ 
Ме{од$, ГАМР, Вазе! — ХийсН, ТМ, Е-2, 1953. 

7. Тротзоп И/ Тгеа4зе оп Мага! РНуозорНу / \/.ТПотбоп апа Р.Та!. 
Раге 1. СатЬчаде Упмег<Ку Ргез$, 1879. 

8. Меркин Д.Р. Введение в теорию устойчивости движения 
/ Д.Р.Меркин. - М.: Наука, 1971. 

9. Бутенин Н.В. Введение в теорию нелинейных колебаний 
/ Н.В.Бутенин, Ю.И.Неймарк, Н.А.Фуфаев. - М.: Наука, 1987. 


Материал поступил в редакцию 18.12.07. 


2АКОМОКОТМУ \.1.’, МСУУЕМ РОМС АМН”*, РНАМ ОТМН ТОМС' 
(Роп За Тесптка! Упмег®Ку* , Козюу- Боп, Кизча; 
пище оЁ Меспапс$"”, Напо!, Мепат) 


ТНЕ ЕМОГОТТОМ ТКАЗЕСТОКУ ЭТАВТЕТТУ РОВ МЕСНАМТСАЁЕ 
ЗУЗТЕМ 1ТМТЕКАСТТМС МЛТН ТВТВОМЕОТА 


183 


Раздел “Машиностроение ” 








Сопаоп$ ог тоЧоп $аБИКу 105$ 15 сопуейпа геаЧ\е фо теспапгса! сует, 
иегасуипа ми тета, Гогттд шт ИБосоирИпа. Сопаюп$ оЁ Гогттд апа гое 
ог 45 $раН\е, дугозсорс, роепНа! ап сисШаг Гогсез, м СН Гогте4 пажигаЙу 
т БосоирИпд апа спапдеа т емоийоп \гапзГогтаНоп, аге зом. 


ЗАКОВОРОТНЫЙ Вилор Лаврентьевич (р. 1940), доктор технических 
наук (1981), профессор (1983), заведующий кафедрой «Автоматизация 
производственных процессов» ДГГУ (1981), проректор по научно-исследо- 
вательской работе ДГТУ. Окончил (1962) РИСХМ по специальности «Авто- 
матические, телемеханические и электроизмерительные приборы и устрой- 
ства» Работает в области динамического мониторинга, взаимодействия 
сложных управляемых систем со средами. 

В 1984 г. совместно с сотрудниками Киевского политехнического института 
и Южмаша стал лауреатом Государственной премии в области науки и тех- 
ники Украины. 

Имеет более 250 научных работ. Руководит аспирантурой и докторантурой. 


НГУЕН ДОНГ АНЬ (р. 1957), доктор физико-математических наук, профес- 
сор, заведующий отделом Института механики академии наук и технологий 
Вьетнама. Председатель Совета Вьетнама по механике. Ученик академика 
АН СССР Ю.А. Митропольского. Область научных интересов: асимптотиче- 
ские методы в нелинейной динамике, статистическая динамика нелинейных 
механических систем. 


ФАМ ДИНЬ ТУНГ (р. 1980), аспирант Донского государственного техниче- 


ского университета. Область научных интересов: эволюция механической 
системы, взаимодействующей с процессом резания. 


184 


